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0> ■ Abstract 

Ch I The position and movement of poles of the amphtude for elastic rj~meson 

^ ' scattering off the light nuclei ^H, ^H, ^He, and ^He are studied. It is found 

■ that, within the existing uncertainties for the elementary -qN interaction, all 

CN , these nuclei can support a quasi-bound state. The values of the r/-nucleus 
scattering lengths corresponding to the critical ryA'^-interaction which pro- 

' duces a quasi-bound state are given. 

in ■ 
cn 

O ■ PACS numbers: 25.80.-e, 21.45. +v, 25.10.+S 

o 

^ I Since meson factories cannot produce r^-meson beams, these particles are available for ex- 
perimental investigations only as products of certain nuclear reactions where they appear as 
final-state particles. Therefore, final-state interaction effects are the only source of infor- 
^ . mation about the ?7-meson interaction with nucleons. In this connection, ?7-nucleus systems 
P5 ! can play an important role in investigating the ryA^-dynamics, especially if they can form 
quasi-bound states. In this case, the final-state 77-mesons can be trapped for a relatively 
long time, and thus the properties of the ?7A^-interaction can be studied. 



X 



Estimations, obtained in the framework of the optical model approach [0,0 , put a lower 
bound on the atomic number A for which an ?7-nucleus bound state could exist, namely 
A > 12. In Ref. [0, the formation of r^-nucleus states has been investigated, using the 
standard Green's function method of many-body problems. There it was found that an 
rj^^O bound state should be possible. Experimentally the cross sections of pion-collisions 
with lithium, carbon, oxygen, and aluminum, however, gave no evidence for the existence 
of T] bound states with these nuclei Q . 

A new theoretical analysis of the problem predicted a binding of the ?7-meson to ^^C 
and heavier nuclei, however, with rather large widths. The formation of an rj'^He bound 
state was studied in a more recent work by Wycech et al. using a modified multiple scat- 
tering theory. These authors obtained a comparatively large negative value for the real part 
of the 77-nucleus scattering length, which was interpreted as an indication that an //-nucleus 
bound state could exist. We note that previous results of ours, concerning the rj scattering 
lengths with light nuclei [0-0, showed that the rj — '^He scattering length can have an even 
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larger (negative) real part than that of Ref . . 



In Ref. a preliminary investigation on the possibility of r^-meson binding in the d, 
t, ^He, and ^ife systems was made within the framework of the Finite-Rank Approxima- 
tion (FRA) of the nuclear Hamiltonian [jlT],|12|. The FRA approach treats the motion of 



the projectile (?7-meson) and of the nucleons inside the nucleus separately. As a result the 
internal dynamics of the nucleus enters the theory only via the nuclear wave function. In 
10| , these wave functions were approximated by simple Gaussian forms, which reproduce 



the nuclear sizes only. In the present work, we perform calculations with more realistic 
nuclear wave functions, obtained via the so-called Integro-Differential Equation Approach 
(IDEA) ||T3Hl7|. We study, in particular, the position and movement of poles of the elastic 
amplitude of r^-meson scattering off the light nuclei ^H, ^H, ^He, and ^ife. 



The approximate few-body equations in the FRA approach enable us to calculate the 
77-nucleus T-matrix 

T{k',k;z) =<k',^o\T{z)\k,tPo> , (1) 

at any complex energy. That is, we can locate the poles of the T-matrix in the complex 
momentum plane p = ^Jl^iz. Here, k is the r^-nucleus momentum, z the total energy of the 
system, /i the ?7-nucleus reduced mass, and '?/'o the nuclear ground-state wave function. 



For the low energies and the light nuclei with only one bound state, being considered, it 
appears justified to approximate the target Hamiltonian Ha by its discrete spectrum 

Ea ~ ^ol^o >< V'ol ■ (2) 

Here iV'o > stands for the ^if, ^if, ^ife, "t/e bound states, respectively, and £^0 for the corre- 
sponding binding energies. 

As a result, we obtain ||^ for the T-matrix the following equation 

nz) = E^°W + E^°(^)l^o > 7 ^ < HTiz) , (3) 

i=i i=i [z - no)[z - Ho - to) 

where Ho is the //-nucleus kinetic energy operator, A the number of nucleons. The T^{z) 
are Faddeev-type components of an auxiliary T-operator, which obey the system of coupled 
equations 

T^{z)=U{z)+U{z)-^-—J:T^{z) . (4) 
{z - ho) 

Here ti describes the scattering of the r^-meson off a nucleon at point r^, where is the 
vector from the nuclear center of mass, which can be expressed in terms of the relative 
Jacobi vectors {r} of the nucleons. In mixed representation, the operator ti is given by 

ti{k', k; r; z) = tr,N{k\ k; z) exp i{k — k') ■ ri 
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with trfN^k', k; z) being the off-shell r]N amplitude. 

Thus, to calculate the T-matrix (|1]) for any fixed value of the complex parameter 
we have to proceed in three steps. First, the coupled integral equations 



Tl'{k',k;f;z) =ti{k',k;f;z) + 



d?k" ti{k', k"; r; z) \ ^ rpQ/j^n r 



(2vr)= 



r; z) 



z — 



(5) 



2/x jj^i 



are to be solved for a number of points r in configuration space, sufficient to perform in a 
second step the integration 



< 



k'Mj:T^{z)%^o >= I d'^^^'^r\Mrl\'j:T^{k\k;f;z). 
i=i ■' i=i 



(6) 



Having determined these matrix elements, it remains, as a final step, to solve the integral 
equation 



T{k',k;z) =<^',^olE^°WI^'^o> +So 



i=l 



^3j^n <k',iJo\ET^{z)\k",^0> 

(2^ (.-f)(.-^o-& 



T{k", k; z) 



2^l 



(7) 



Note that after partial-wave decomposition both equations and (0) become one- 
dimensional. As an input information, we need the ground-state wave functions ipo of the 
nuclei involved and the two-body T-matrix t^N- 

The "00 (r) ^=3 and A = 4 were obtained by means of the IDEA [|T^]. In this method 
the A-body bound state wave function is expanded in Faddeev-type components. 



given as solutions of 



i<j<A 



k<l<A 



where fij = fi — fj. The IDEA is then introduced using the ansatz 

V^,,(f) = i^[,„](f)P(o„p)/p(^-^)/^ 



(8) 



(9) 



(10) 
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with p = 2/AJ2''"ij being the hyperradius, D = 3{A — 1), and H[L^]{r) the harmonic 



polynomial of minimal degree [L^] [T^. For [L^] = the IDEA reads 



P 



k<l<A 



(11) 
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where Vo(p) is the so-called hypercentral potential [|T8[. Projecting Eq. (0) onto the 



space provides us, for spin dependent nucleon-nucleon potentials, with two coupled inte- 
gro differential equations for the symmetric S and mixed symmetric S' components of the 
function P^{Qj, p) , n = S , S' . More details and explicit equations are found in Refs. |13|,|15 



For the nuclear ground states, we use the fully symmetric S-wave components obtained 
with the semi-realistic Malfliet-Tjon I-III (MT I-III) nucleon-nucleon potential |T^. The 
corresponding two-, three-, and four-body binding energies are 2.272 MeV, 8.936 MeV, and 
30.947 MeV, while the root mean square (r.m.s.) radii are 1.976 fm, 1.685 fm and 1.431 
fm, respectively. The omission of Coulomb effects and of the mixed symmetry components 
makes and ^He indistinguishable. In order to compensate partly for this omission, we 
use in Eq. (|^) the experimental values for masses and binding energies of the nuclei ||20|| . 



At low energies the t]N interaction is dominated by the iV*(1535) Sn - resonance. For 
the T^A^-amplitude we, therefore, choose the separable form 



A 



+ a^){z -Eo + iT/2){k^ + a2) 



(12) 



with Eq = 1535 MeV - {rriN + nir,) and F = 150 MeV 



a, we use the results of Refs. p2|j23[| . There the same rjH 
was employed, and a was determined via a two-channel fit to the txN 
experimental data. 



. To find the range parameter 
N* vertex function [k"^ + a^^-i 
ttN and vriV 77N 



Three different values for the range parameter a are available in the literature, namely, 
a = 2.357 fm-i 1221, a = 3.316 fm'^ 121, and a = 7.617 fm"^ l2l. Since there is no 



^ ^ a = 3.316 fm~ ||2^, and a ^ ^ 

criterium for singling out one of them, we use all three in our calculations. The remaining 
parameter A is chosen to provide the correct zero-energy on-shell limit, i.e., to reproduce the 
TjN scattering length o^at. 



t^^(0,0,0) 



27r 

PriN 



'"qN- 



(13) 



Different analyses provided values for the real part Rear,N in the range 0.27 0.98 fm and 
for the imaginary part Imar^N in the range 0.19-i-0.37fm 1^^. To examine at which value of 
a^TV within the above ranges an 77-nucleus bound state exists, we parametrize the scattering 
length as follows 



(^riN = ((70.55 + ig'0.30) fm, 
where g and g' are adjustable parameters. 



(14) 



Since arjN is complex, the ?7-nucleus Hamiltonian is non-Hermitian and its eigenvalues are 
generally complex. In this case, eigenvalues attributed to resonances and quasi-bound states 
are located in the second-quadrant of the complex p-plane |]2B|. The energy Eq = pl/2fi 
corresponding to a pole at p = po, 



En = 

° 2/i 



(RepoY - {Impof + 2i{Repo){Impo] 



(15) 
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has a negative real part, RcEq < 0, only if pq is above the diagonal of this quadrant. Such 
a pole is related to a quasi-bound state. For pq below the diagonal we have Re Eq > and 
the pole is then attributed to a resonance. Therefore this diagonal is critical: when crossing 
it from below a resonance pole becomes a quasi-bound pole. 

Fixing g and g' of Eq. ([I4D to g = g' = 1 and varying the complex momentum p = -^/S/ii^, 
we located the poles close to the origin p = 0. The results obtained are given in Table 1. For 
one choice of the range parameter, namely a = 2.357 fm~^, the positions of the poles found 
are shown by the open circles in Fig.l. It is seen that for the rjd, rjt, and rj^He systems, these 
poles lie in the resonance region, while for the rj'^He system the pole is in the quasi-bound 
region. 

Increasing g while keeping g' = 1, the resonance poles are moving towards, and finally 
cross, the diagonal. In the deuteron case, the corresponding trajectory is depicted in Fig. 1 
by the solid curve which crosses the diagonal when g = 1.6536. 

To find the relationship of poles above the diagonal to physical bound states, we grad- 
ually removed the imaginary part of a^^v by fixing g and decreasing g' in Eq. (^) to zero. 
The imaginary part of the Breit-Wigner factor in Eq. (0) was also decreased, using the 
same parameter g', so that it goes over into {z — Eq + ig'T /2)~^. For g' = Q the Hamiltonian 
becomes Hermitian and, hence, the bound state poles in this case must be on the positive 
imaginary axis. The dashed curve in Fig. 1 is the trajectory of the r}d bound state pole 
(with g = 2) when g' decreases from 1 to 0. It is seen that the final position of the pole lies 
on the positive imaginary axis. This supports our interpretation of poles above the diagonal 
as quasi-bound states. 

By varying the enhancing factor g for each of the 77-nucleus systems under consideration, 
we found the values which generate quasi-bound states on the diagonal. They are given in 
the Table 2. These values correspond to an rjN attraction, which generates an //-nucleus 
binding with RcEq = 0. Further increase of g moves the poles up and to the right, enhanc- 
ing the binding and reducing the widths of the states. The value of Re a^^r which provides 
the critical binding lies within the range [0.27, 0.98]fm used in the literature. Therefore, an 
?7-nucleus quasi-bound states may exist for A > 2. If this is not the case, then at least 
near-threshold resonances (poles just below the diagonal) could exist. However, as can be 
seen in Table 1 and 2, the widths of such quasi-bound and resonance states could be small 
only for the rj'^He system. 

In Table 3 we present the //-nucleus scattering lengths calculated with parameters gen- 
erating the critical ?7-nucleus binding. From this table we see that the real part of the 
7/-nucleus scattering length can be small despite the existence of a quasi-bound state. This 
is due to the non-Hermitian nature of the rjN interaction. Being complex, this interaction 
generates critical poles rather far from the origin and their influence on the scattering length 
(the value of the amplitude at the origin) is not very strong. 

In conclusion, we have shown that the uncertainties in the riN scattering length allow 
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for choices of parameters in the ryA^-amplitudc that may generate poles in the 77-nucleus 
amphtudes considered, which can be attributed to quasi-bound states. 
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Table 1. Positions of poles po — y/2jlEQ of the 77-nucleus amplitudes with g — g' — 1 ior 
the three values of the range parameter a. 





Po (fm ^) 


Eq (MeV) 


a (fm-i) 




—0 90259 + 70 35870 


'.} 1. . t: '-.J \J v^U>\JU±. 


2 357 


T] d 


-0.8 159 1 + /0. 32195 


28.061 - /21.976 


3.316 




-0.82460 + i0.30423 


26.935 - 723.006 


7.617 




—0 56045 + iO 23859 


10 906 _ i\\ 341 


2.357 


Tjt 


-0.55511 + 70.26826 


10.015 -712.630 


3.316 




-0.51725 + i0.27896 


8.0456 -^12.238 


7.617 




-0.54692 + 70.24478 


10.143 -7II.354 


2.357 


rj^He 


-0.50815 + 70.30402 


7.0305 - 713.102 


3.316 




-0. 18310 + /0. 339 18 


5.0099 - /■13.909 


7.617 




-0.16504 + i0.27876 


-2.0540 - 73.7447 


2.357 


r) ^He 


-0.20215 + i0.38726 


-4.4403 - i6.3718 


3.316 




- 0.25931 + 70.45846 


-5.8175-79.6766 


7.617 



Table 2. The parameter g generating the jy-nucleus amplitude poles po — -\/2/lB^ on the 
diagonal for the three values of the range parameter a and g' — 
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Po (fm ^) 


Eo (MeV) 


a (fm-^) 


77 d 


1.6536 


-0.32527 + 70.32527 


-79.7026 


2.357 


1.5605 


-0.33541 + 70.33541 


-710.317 


3.316 


1.5260 


-0.33670 + 70.33670 


-7IO.397 


7.617 


Tjt 


1.3624 


-0.33515 + ^0.33515 


-i9.5266 


2.357 


1.3055 


-0.35190 + i0.35190 


-ilO.503 


3.316 


1.2436 


-0.35186 + 70.35186 


-710.500 


7.617 


rj^He 


1.3306 


-0.34034 + 70.34034 


-79.8239 


2.357 


1.2171 


-0.36267 + 70.36267 


-711.155 


3.316 


1.1421 


-0.37631 + 70.37631 


-7I2.OIO 


7.617 




0.86222 


-0.20641 +i0.20641 


-i3.4679 


2.357 


0.80813 


-0.26522 + z0. 26522 


-i5.7255 


3.316 


0.79578 


- 0.35215 + 70.35215 


-710.094 


7.617 



Table 3. The T^-nuclcus scattering lengths for the parameter g of Table 2, which generate 
the condition for binding (BeE — 0). 





a = 2.357 (fm^^) 


a = 3.316 (fnr^) 


a = 7.617 (fm^^) 


Tj d 


0.171 + 75.99 


-0.198 + 74.57 


-0.318 + 73.52 


r)t 


-3.65 + 73.49 


-2.91 + 73.02 


-2.19 + 2.70 




-3.49 + i3.67 


-2.66 + i3.31 


-1.96 + i2.86 




-3.43 + 72.60 


-2.81 + 72.14 


-2.30 + 71.72 
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FIGURES 




FIG. 1. The ry-nucleus elastic scattering amplitude pole positions in the complex p-plane. 
The open circles correspond to g=l. The solid curve is the r/d-amplitude pole trajectory when g 
increases from g=l to g=2. The dashed curve shows the trajectory of the r)d pole with g=2 and 
with g' varied until the r)N interaction becomes real. 
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